Family of analytic entanglement monotones 
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We derive a family of entanglement monotones, one member of which turns out to be the negativ- 
ity. Two others are shown to be lower bounds on the I-concurrence, and on the I-tangle, respectively 
[P. Rungta and C. M. Caves, to appear in Phys. Rev. A]. We compare these bounds with the I- 
concurrence and I-tangle on the isotropic states, and on rank-two density operators resulting from 
a Tavis-Cummings interaction. Our results provide a global structure relating several different en- 
tanglement measures. Additionally, they possess analytic forms which are easily evaluated in the 
most general cases. 
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Entanglement has proven to be a useful resource in 
the implementation of quantum information processing 
protocols jjj. This observation, coupled with the impor- 
tant role played by entanglement in the foundations of 
quantum theory, has led to the search for a quantitative 
theory of entanglement. In this context, a deep connec- 
tion between the concepts of separability and of positive 
maps has been established ||, and several different 
measures of entanglement have been proposed Q . 

The concurrence and the tangle are two related en- 
tanglement measures yielding analytic forms in the case 
of two qubit quantum systems [|| . Due to their intimate 
connection with the entanglement of formation, and with 
the phenomenon of entanglement-sharing ||, they have 
proven to be particularly useful tools for studying funda- 
mental issues in quantum mechanics. Recently, analytic 
forms for these entanglement measures have been found 
for bipartite quantum systems of arbitrary dimensions in 
an overall pure state . However, the generalizations of 
these quantities, known respectively as the I-concurrencc 
and the I-tangle ^ , to mixed states of a bipartite system 
of arbitrary dimensions, involves a difficult minimization 
over ensemble decompositions. So far, it has been pos- 
sible to calculate analytically the I-concurrence and the 
I-tangle only for the isotropic states ||, and in the spe- 
cial case that the density matrix of an arbitrary bipartite 
quantum system has a rank no greater than two Q . 

Except for the cases mentioned above, analytic forms 
for the I-concurrence and the I-tangle are lacking. In 
this work we present a new family of entanglement mono- 
tones, certain members of which constitute easily com- 
putable lower bounds for the I-concurrence and for the 
I-tangle. In order to construct these functions we make 
use of the fact that the I-concurrence and the I-tangle are 
convex-roof extensions of the generalized concurrence for 
pure states of arbitrary dimensions, and of the square of 
this quantity, respectively Therefore, the I-concurrence 
(I-tangle) can be bounded from below by any convex 
function which agrees with the I-concurrence (I-tangle) 
on the set of bipartite pure states. This requirement 
is seen to be satisfied by specific members of the set of 



entanglement measures introduced herein. Each is a sim- 
ple function of the negative eigenvalues generated via the 
partial transposition operation, and may be easily calcu- 
lated with any standard linear algebra package. 

In order to estimate the quality of these functions as 
lower bounds, we compare their values with the values 
of the I-concurrence and the I-tangle on the family of 
isotropic states. The I-tangle, and the corresponding 
lower bound, are also compared numerically for rank- 
two density matrices arising in the context of the Tavis- 
Cummings model fto| . 

The construction that follows is based on the theory 
of majorization. This formalism has become a very use- 
ful tool in characterizing the relationships among den- 
sity matrices, ensemble decompositions, and measure- 
ment processes [ pd| , and has led to new insights in the 
structure of quantum algorithms [T^| , and in the problem 
of Hamiltonian simulation |k| . 

The following is a brief review of the main tenets of 
majorization theory. The reader is referred to Jl4| for 
extensive background on the subject. Given two vectors 
x and y in R™, we say that the vector x is majorized by 
the vector y, denoted by the expression x -< y, when the 
following two conditions hold: 
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Here, the symbol J, indicates that the vector coefficients 
are arranged in decreasing order. 

Majorization is naturally connected with the idea of 
comparative disorder p^ j. In fact, x -< y if and only 
if there exists a doubly stochastic matrix A such that 
x = Ay. A matrix A is doubly stochastic if its coefficients 



are non-negative and 



o>ik 



Yji a kj = l,Vfe. If we 



consider x and y to be probability distributions, then the 
fact that x is majorized by y expresses the idea that x is 
more disordered, in a quantifiable sense, than y. 
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In the case that only condition (|l|) holds, we say that 
x is weakly submajorized by y. This is denoted by the 
expression x -< w y. We will make use of the following 
two results concerning majorization [ p"4| : 

x -< w y e R" =*> x + -< w y + (3) 
x^ w yeK'l => xP^ w yP,Vp > 1, (4) 

where the operations x p and x + act on each component 
of x individually. The x + operation simply converts each 
of the negative entries in x into a zero. 

The following relation allows us to construct a family 
of convex functions of the negative eigenvalues of an her- 
mitian matrix. Given any two hermitian matrices A and 



\{A + B) -< X l (A) + X l {B), 



(5) 



where A(T) denotes the vector whose coefficients are 
eigenvalues of T. Let us now define the vectors X(T) = 
— A(T) = A(— T), such that the negative coefficients in 
A(T) become positive in A(T). Clearly Eq. (||) also holds 
for the vectors A(T), i.e., 



\{A + B) -< X l (A) + X l (B). 



(0) 



Recognizing that the coefficients of X(T) belong to R", 
and applying property (|J), Eq. (||) becomes 

(X(A + B)Y^ W (\^A) + ~X^B)Y . (7) 

The coefficients of the vectors in Eq. (0) are, by defini- 
tion, members of K™ . Thus, using property (0) we obtain 



(x(A 



+ B] 



(x l (A) + X l { 
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According to condition (|l|) for k = n, we have 



n n r 

]T(A+L4 + B)) P <£ (\l(A) + \j(B) 
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i=l 



(8) 



, (9) 



where we have removed the ordering of the vector on the 
left hand side. The term inside square brackets on the 
right hand side of Eq. (|^) can be bounded from above 
by X L +(A) + X l +(B), yielding 



n n 

J2 (A+(A + B)) P < (Al + (A) + Aj+(B) 



(10) 



Taking the p-th root of Eq. (^0|), and using Minkowski's 
inequality p4[ , we obtain 
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(11) 



The terms in square brackets on the right hand side of 
Eq. ( |TT| ) are the sums of the positive coefficients of A(^4) 
(A(-B)) to the p-th power, or equivalently, to the sums of 
the absolute values of the negative coefficients of X(A) 
(A(-B)) to the p-th power. Thus, we see that the quantity 



M, 



p(a) = ( Yl \ X ( A )\ P 



i/p 



(12) 



A(,4)<0 



obeys the triangle inequality on the set of hermitian ma- 
trices for any p > 1. In particular, M p (A) is a convex 
function, i.e., 

M p (aA + j3B) < aMp(A) + (3M P (B) (13) 

for a and (3 in the interval [0, 1] such that a + (3 = 1. 
Similar results may also be shown to hold for the set of 
functions 



N p (A) = M p (A) p y P > 1. 



(14) 



The I-concurrence C{p) of a density matrix p acting 
on a bipartite d-dimensional Hilbert space Ti. is defined 

by | 

C(p) = { min } j£>C(*0|p = l>l**><**l} ' ( 15 ) 

with the concurrence C(^) of a bipartite pure state 
given by 



2 2 



(16) 



The Ci's in Eq. (16) are the coefficients of the state 
\^>), written in the Schmidt decomposition Q|. The I- 
concurrence is the convex-roof extension of C(JS), and 
represents the average value of the pure state concur- 
rences for an ensemble decomposition of p, minimized 
over all possible ensemble decompositions. The I-tangle 
has a similar construction, and is given by 



r(p) = min 



(17) 

Due to a result by Uhlmann |l5| , the I-concurrence and 
the I-tangle are known to be the largest convex functions 
defined on the set of density operators which agree with 
C(^) and C 2 (\I/), respectively, over the set of bipartite 
pure states. Therefore, if we are able to find a convex 
function which agrees with C('J') (C 2 ( x f')) over the set of 
bipartite pure states, then it will automatically constitute 
a lower bound for the I-concurrence (I-tangle). This can 
be accomplished by observing that the coefficients in Eq. 
( p6| ) are the absolute values of the negative eigenvalues 
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of the partial transpose of Since the partial trans- 
position operation is linear, the function M^(p) defined 
by 



Ml\p) = 2M 2 (pP t ) 



(18) 



is a convex function which agrees with the I-concurrence 
on the set of bipartite pure states. Consequently, \p) 
is a lower bound for the I-concurrence, i.e., 



Mf(p) < C(p). 



(19) 



In a similar way it can be shown that the function 
N^ip) = [Mf (p)] 2 is a lower bound for the I-tangle, 



i.e., 



N?(p) < r{p). 



(20) 



These results hold for arbitrary dimensional bipartite 
quantum systems. 

A similar result was obtained in Jl6| for the case of two 
qubits where d = 4. Specifically, it was shown that the 
negativity Af(p) is a lower bound on the concurrence [fi"7| . 
The negativity is defined to be the absolute value of the 
sum of the negative eigenvalues of the partial transpose 
of p. Thus, the negativity is seen to be one member of our 
family of entanglement monotones, i.e., N(p) = Mf (p). 
Additionally, in the case d = 4 the partial transpose 
of p has at most one negative eigenvalue, implying that 
M^^p) also reduces to the negativity in this situation. 

The lower bounds given by Eqs. @ and (§(]) are func- 
tions of the negative eigenvalues produced via the partial 
transposition operation. Hence, they are entanglement 
monotones in their own right. This follows directly from 
the monotonicity of Af(p) fllffl. Specifically, it may be 
shown that monotonicity is preserved in Eqs. ( p^ ) and 
( |l4| ) for choices of p other than one. Further, these quan- 
tities have the additional advantage that they may be 
evaluated in a straightforward manner with the help of a 
standard linear algebra package. 

It has been shown that the positive partial transposi- 
tion criterion is a necessary and sufficient condition for 
separability for d < 6 [|J. In higher dimensions, positivity 
under partial transposition is a necessary, but not suffi- 
cient, condition for separability. However, this is not a 
serious drawback for the usefulness of these lower bounds. 
Indeed, it has been shown by numerical experiments and 
theoretical results that the volume of the set of density 
operators with positive partial transpose decreases expo- 
nentially with the dimension d of the Hilbert space jl8) . 

The exact values of the I-tangle for the isotropic states 
pp was analytically calculated in ||. Isotropic states 
describe a quantum system composed of two subsystems 
of equal dimension d. They are mixtures formed by the 
convex combination of a maximally mixed state and a 
maximally entangled pure state, i.e., 



PF = (1 



(21) 
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FIG. 1: Comparison of r (pf) and (pf) as functions of 
the fidelity F for different dimensions d. Solid line: r (pf) for 
d = 100. Dashed line: Af ' (p F ) for d = 100. Dot-dashed line: 
t (p F ) for d = 3. Dotted line: iVf ' (p F ) for d = 3. 



Here, / is the identity operator acting on a d-dimensional 
Hilbert space, and is the state given by 



1 



(22) 



The parameter A in Eq. (|Tj) can be related to the fidelity 
F of pp with respect to the state (where F = 

(*+|/9 F |*+) e [0, 1]), via the relation 



A = 



d 2 F- 1 
d?-l ' 



(23) 



It has been shown that the isotropic states are separable 
for F < 1/d |9j. 

The value of the lower bounds M^{p) and A/J (p) on 
the isotropic states can be calculated easily. Since the 
partial transposition operation is linear, and since the 
identity operator is invariant under this operation, the 
eigenvalues of p P p are readily obtained. They are given by 
(1 — X)/d 2 ±\/d with multiplicity d(d±l)/2, respectively. 
The negative eigenvalues (1 — A) jd 2 — X/d become positive 
when A < l/(rf+l), or equivalently, whenF < 1/d. Thus, 



Mf(p F ) 



(A_i +A) ,yfciT A >i/(d+i) 

A<l/(d+l) 



and 



N?( PF )=[M$ t (p F )Y 



(24) 



(25) 



The behaviors of the I-tangle r(p F ) and of N^^Pf) 
for the isotropic states are depicted in Fig. [I]. For 
d = 2, the two functions assume the same values, while 
for larger dimensions and constant fidelity, the difference 
between the lower bound and the I-tangle increases. In 
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FIG. 2: Comparison of r (p a /) and iVj (p a /) calculated for 
the two-atom TCM as functions of the effective time <?i. The 
upper curve shows the evolution of r(p a f), and the lower 
curve demonstrates the similar qualitative behavior of the 
lower bound N% 1 (paf)- 



the limit d — ► oo, t(pf) and N$ (Pf) behave as \/2F and 
2F 2 , respectively. Similarly, the I-concurrence C(pf) 
and Mf'(p F ) may be calculated analytically. Here we 
find that the two quantities agree over the isotropic states 
for any dimension d, demonstrating that this is a tight 
lower bound. 

It has been shown |l for the case of rank-two density 
matrices that an analytic formula for the I-tangle ex- 
ists. This class of density matrices arises naturally in the 
context of the two-atom Tavis-Cummings model (TCM). 
The TCM describes the interaction of an ensemble of N 
two-level atoms with one mode of the quantized electro- 
magnetic field in the dipole and rotating wave approxi- 
mations Jk}. In the case where N — 2 and the initial 
state of the overall system is pure, tracing over one of 
the atoms results in a reduced density operator for the 
remaining atom and field subsystem of rank at most two. 
Consequently, the entanglement between the field and 
the atom can be quantified using Osborne's formula for 
the I-tangle §. 

Fig. H shows the evolution of r (p a f), and of N^ipaf), 
as functions of the effective time gt, where g represents 
the coupling strength between the two-level atoms and 
the electromagnetic field, and p Q / is the reduced density 
operator for one of the atoms and the field. We have 
considered an ensemble of two atoms, both of which are 
initially in the excited state. The field is initially in a 
coherent state with an average photon number (n) = 100. 
The I-tangle and its lower bound show clear differences 
in their magnitudes. However, the lower bound preserves 
the basic structure present in the evolution of the tangle. 
This agreement is a direct result of the monotonicity of 

The entanglement measures analyzed above belong to 
larger classes of monotones defined by taking A — p pt in 



Eqs. (|12j) and (J14|). The negativity defined for a system 
of two qubits is seen to be one member of this set. Other 
instances correspond to lower bounds on the quantities 
I-concurrence and I-tangle, which are useful tools for in- 
vestigations of quantum information theoretic concepts, 
and of fundamental quantum mechanics. Apart from of- 
fering a larger structure from which to view these differ- 
ent entanglement measures, this new family of functions 
also possesses analytic forms which are easily computed 
even in the most general cases. It will be interesting to 
see if further connections between seemingly unrelated 
measures of entanglement can be found using this for- 
malism. 
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